


On-chip low loss heralded source of pure 
single photons 

Justin B. Spring! *, patrfck S. Salter^, Benjamin J. Metcalf^ Peter C. 

Humphreys^ Merritt Moore ^ Nicholas Thomas-Peter', Marco 

Barbieri', Xian-Min Jin'^, Nathan K. Langford' ''j W. Steven 

Kolthammer', Martin J. Booth^^, and Ian A. Walmsley' * 

^Clarendon Laboratory, University of Oxford, Parks Road, Oxford 0X1 3PU, United Kingdom 
^f^ ^Department of Engineering Science, University of Oxford, Parks Road, 

^-H Oxford 0X1 3PJ, United Kingdom 

f^ ^ Department of Physics, Shanghai Jiao Tong University 

^sj Shanghai 200240, People's Republic of China 

^ Department of Physics, Royal Holloway, University of London, 
London TW20 OEX, United Kingdom 
. ' Centre for Neural Circuits and Behaviour, University of Oxford, Mansfield Road, 

<^ Oxford 0X1 3SR, United Kingdom 

*Corresponding authors: j. spring 1 @ physics. ox.ac. uk, i. walmsley 1 @ physics. ox.ac. uk 

(N 

^^ Abstract: A key obstacle to the experimental realization of many photonic 

C~| quantum-enhanced technologies is the lack of low-loss sources of single 

Qh photons in pure quantum states. We demonstrate a promising solution: 

generation of heralded single photons in a silica photonic chip by spon- 
taneous four-wave mixing. A heralding efficiency of 40%, corresponding 
to a preparation efficiency of 80% accounting for detector performance, 
^H is achieved due to efficient coupling of the low-loss source to optical 

'— ' fibers. A single photon purity of 0.86 is measured from the source number 

. statistics without filtering, and confirmed by direct measurement of the joint 

^ spectral intensity. We calculate that similar high-heralded-purity output can 

1— I be obtained from visible to telecom spectral regions using this approach. 

OO On-chip silica sources can have immediate application in a wide range of 

r single-photon quantum optics applications which employ silica photonics. 
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1. Introduction 

Photonic quantum-enhanced technologies aim to employ nonclassical states of light to surpass 
classical performance limits in diverse fields including computation |[Tl|2l, metrology ||3|, and 
communication |4|. An impediment to faster progress is the quality of available single photon 
sources. Building the low-loss sources of high purity single photons necessary for quantum- 
enhanced performance has proven challenging 15] |6] |7l |8] HI . 

Heralding spontaneous emission from a nonlinear-optical material has to date been the most 
common method of generating single photons |fT0l[m[T2l|T3J . Nonlinear processes such as 
spontaneous parametric down-conversion (SPDC) or spontaneous four-wave mixing (SFWM) 
can be used to create pairs of photons. Due to this pairwise emission, detection of one photon, 
the heralding photon, indicates the creation of its partner, the heralded single photon. A key 
source metric is the preparation efficiency rjp, the conditional probability that a heralded photon 
is delivered to its application given detection of the heralding photon. For example, the rate at 
which multiphoton states are generated from single-photon sources scales exponentially with 
rjp, even with multiplexing strategies which ideally achieve near-deterministic emission lfT4ll . 
In numerous applications, rjp is a crucial parameter for a quantum method to demonstrate true 
advantage over a classical approach |l6]|2l|8]|9l. 

Four-wave mixing in a silica waveguide is a promising route to achieving exceptionally high 
rjp lfT2l [TSJI . In heralded photon sources, reduction in rjp results from loss of the heralded 



photon, due to scattering or imperfect mode-matching at interfaces. Waveguides in commonly 
available silica glasses minimize these effects due to their exceedingly low optical loss and ex- 
cellent mode-matching to optical fiber, a ubiquitous component in quantum photonics. While 
silica's relatively small x nonlinearity affects the emitted photon flux, in many applications 
it is loss that is fundamental to demonstrating quantum enhancement, not flux. In fact, for her- 
alded single photon sources one must deliberately keep emission rates low to avoid unwanted 
heralding of more than one photon. The transverse confinement of the waveguide allows these 
desired photon production rates to be reached at readily available pump powers. 

Typical quantum applications require heralded photons in pure quantum states in addition 
to high source rjp. In general, however, the heralded pair source generates mixed quantum 
states. Energy and momentum conservation can lead to entanglement between multiple spatial 
and spectral modes of the emitted photon pairs fT6^, TT"]. If the heralding photon is detected 
but its mode is not resolved, then the heralded photon is left in a mixed quantum state of all 
possible modes. One approach to achieve a high-purity heralded photon is to use spectral and 
spatial filters on the heralding field which ensures the detector responds to only a single spatio- 
temporal mode. Such filters reduce the rate at which heralded photons are emitted. On-chip 
SPDC sources are capable of high photon flux lITSl IT9l l20l 1211 that allows acceptable count 
rates even when filters are used to achieve good photon purity. For SFWM in silica, on the 
other hand, the reduced count rate from this approach could lead to unacceptably long data 
collection times. 

An alternate approach we employ here is to engineer the source to emit photons into a single 
pair of modes. As a consequence, the heralding and heralded photon are not entangled. The 
resulting factorable output allows heralded photons of high purity without filtering. Guided- 
wave photonics allows the precise control of optical modes needed to construct such a source. 
Previous silica sources have demonstrated this strategy using optical fibers. In these structures, 
however, fabrication inhomogeneity llTni22l or sensitivity to the local environment llT2ll23ll24l 
|251 prevented a robust scalable solution. While on-chip SFWM sources have been demonstrated 
in chalcogenide glass ll26l and silicon [271, these devices exhibited relatively high loss and did 
not pursue the factorable source design strategy described here. 

Here we report the first photon source on a silica photonic chip. Our SFWM source generates 
heralded single photons of purity P = 0.86 and rjp = 80% without any filtering. Heralded single 
photons are emitted at a rate of 3.1 • 10^ photons per second with a pump field power of 150 
mW. We use a birefringent waveguide to phasematch SFWM at frequencies where the red- 
detuned photon lies far beyond silica's Raman gain peak. This minimizes spontaneous Raman 
scattering ||28l[T2l[T5l . which is the principal noise source in many silica sources |29, 30 1. A 
solid-clad silica waveguide provides a convenient and robust platform for our source which we 
foresee finding immediate application in integrated quantum optics experiments that frequently 
rely on similar integrated silica architectures El IT3] [3111321 [33] [34|[35ll . 

2. Experimental overview 

Our source uses a 4 cm long waveguide fabricated by femtosecond laser writing in an undoped 
silica chip (Lithosil Ql) |[36l[37l . We use adaptive optics to shape the writing beam |38l and 
produce an elliptical transverse mode yielding a birefringence of A« = 10^^. By measuring the 
insertion loss and imaging the spatial mode, we determine the propagation loss in our waveg- 
uide to be less than 0.4 dB/cm r38 1. 

A single pump field with central wavelength Xp = 729 nm is generated by an 80 MHz Ti- 
Sapphire oscillator whose spectral bandwidth AXp is adjusted with a mechanical slit in a 4-f 
line, as shown in Fig. [it a). We describe the high (low) frequency photon in the emitted pair as 
the signal (idler). A dichroic beam splitter (Semrock FF740) separates these signal and idler 
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Fig. 1. The experimental layout consists of (a) state preparation and (b-d) characterization, 
(a) A 4-f line with an adjustable slit is used to control the bandwidth of the pump (AAp). 
After the chip, the pump is filtered out while the signal and idler fields are separated by 
a dichroic mirror (DM) and coupled into separate single-mode fibers (SMFs). (b) Source 

count rate and cross-correlation g],- (0) are measured with the SMFs coupled directly into 

(2) 
avalanche photodiodes (APDs). (c) The autocorrelation g^^ (0) is measured by connecting 

the idler fiber to a fiber beam splitter (BS) with a reflectivity of 50% and ignoring the sig- 

nal field, while the signal field is used as a herald when determining g)^ (0). We measure 

f2) 
& (0) by inserting the signal fiber into the BS and ignoring the idler field (not shown), (d) 

To measure joint spectral intensities, the source SMFs are connected to separate monochro- 

mators with multi-mode fibers (MMFs) at the output performing a raster scan over the joint 

spectral range while APDs count in coincidence. 



photons which have central wavelengths of A, = 676 and A,, ~ 790 nm. A combination of spec- 
tral (Semrock 684/24 and 800/12) and polarization filters extinguish the pump field before 
the signal and idler modes are coupled into single-mode fibers. Silicon avalanche photodiodes 
(APDs) (Perkin Elmer SPCM-AQ4C) and an FPGA-based coincidence counter (XiUnx SP-605) 
are used to measure marginal and joint photon statistics as shown in Fig.[TJb)-(d) and discussed 
below. 



3. Nonclassical emission and heralding of single photons 

We first confirm nonclassical operation of our source by measuring the cross- and auto- 
correlations of the signal and idler modes with the setup shown in Fig.[T|b)-(c). Classical fields 

must satisfy the Cauchy-Schwarz inequality (g„- (0))^ < gL (0) ■ g\i (0) where gxy {'^) is the 



second order coherence between modes x and y at relative time delay T 
source, we calculate 

^^"^' ^"^ " A^,A^, 



. For our pulsed 
(1) 



where A^^ {Nxy) correspond to single (temporally coincident) detection events on mode x {x 
and y). The number of trials is given by Np, the number of pulses from the Ti-sapphire pump. 

(2) 

Autocorrelations are described by ^1,^2(0) where x\2 refer to the two output ports of a fiber 
beamsplitter with a reflectivity of 50% as shown in Fig.flTc). 

With the pump bandwidth AA,, = 3.1 nm and 100 mW average power, we measure g\^ (0) = 
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Fig. 2. The count rate for the signal mode alone (Ns) and in temporal coincidence with the 
idler mode (Nsf) is shown as a function of the pump power (blue circles). The heralding 
efficiency, which includes APD detector losses, is calculated according to rjjj = N^i/Ns 
(green). At low pump powers, Tj// decreases due to the heightened importance of detector 

dark counts, which cause false herald events. The heralded autocorrelation of the idler 

(2) 
photon with the signal field as herald, ^^^ (0), increases linearly due to spontaneous Raman 

scattering (red squares). Quadratic (blue) and linear (red) fits to the data are shown with 

(2) 
dotted lines. Error bars are smaller than the markers for the count data and g)^ (0). 



73.5 ± 1.1, g[^\Q) = 1.82 ±0.03, and glf\o) = 1.26 ±0.02 with no background subtraction. 
The Cauchy-Schwarz inequahty is violated here by 49 standard deviations which demonstrates 
a nonclassical correlation between signal and idler modes in the photon number basis. 

It is this correlation in photon number that makes such sources suitable for heralding single 
photons. Using the signal photon as the heralding photon, we measure the heralding efficiency 
iffj =Nsi/Ns- We calculate rjn > 40% for pump powers of 75 to 150 mW as seen in Fig.l2] Since 
rjH is limited by the detector efficiency, which is not our concern here, we estimate rjp, which 
corresponds to the probability that the idler photon arrives at its detector given detection of 
the heralding signal photon. Using the manufacturer's specified detector efficiency, we estimate 
rjp = 80%. The remaining inefficiency is primarily due to loss from coupling to single mode 
fiber, reflection at interfaces, and scattering. Neither the chip nor the fibers were AR-coated. 

We quantify the suppression of undesired higher order photon emission via the heralded 
second-order correlation at zero relative time delay 



8^\0): 
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for which all counts are conditional on detecting a heralding signal photon as shown in Fig. [Tic). 

(2) 12) 

An ideal single-photon source would give g[j (0) — 0. We measure g^ {0) — 0.0092 ±0.0004 
for a pump power of 25 mW. Spontaneous Raman scattering is the primary noise source, which 

(2) (2) 

explains the linear increase in g[j {0) as the source is pumped harder. Our low g)/ (0) values 
are comparable to other SFWM sources ITSl l25l that take advantage of birefringent phase 



matching [28 1, which allows significant suppression of this Raman noise which has inhibited 
other SFWM sources 



4. Controlling the heralded state purity 

The pairwise emission that allows heralding of single photons is not sufficient for many applica- 
tions; these photons must also be in pure states. Undesired correlations between the signal and 
idler fields that are not resolved by the heralding detector instead result in heralding of mixed 
states. Such correlations imply multimode emission in the frequency, polarization, or spatial 
degrees of freedom. Phasematching constraints generally force the signal and idler fields to 
be emitted into single polarizations. Furthermore, our waveguide constrains each of the three 
fields (pump, signal, idler) to a single spatial mode. This is more readily achieved for a SFWM 
source, in contrast to SPDC, due to the reduced disparity in field frequencies. 

We focus on the remaining possibility that correlations are generated in the spectral degree 
of freedom. The effective SFWM Hamiltonian can thus be approximated as lfT6l 

^SFWM == C // dCO,dCOif{CO,, (Oi)a'^ {CO,)P{COi) + H.c. (3) 

where the joint spectral amplitude f{c0s,C0i)= fdco'a{co')a{cOs+cOi—co')(p{cOs,cOi) is a func- 
tion of the pump envelope a{(o) and phasematching function (p{cOs,COi), while ^ depends on 
both the pump intensity and magnitude of the x nonlinearity. We approximate a{co) as a 
Gaussian function and specify AXp as the full-width at half maximum of |a(a))p. The phase- 
matching function results from integrating over the length of the guide (Oi , 0),) — Jq e'^^dz °= 
e'^^sinc(AA:L/2), where the wavevector mismatch is Ak — 2kp — kg — A;,-. We neglect the phase 
mismatch arising from the pump pulse intensity, since this is small for our source. In the normal 
dispersion regime, SFWM is phase-matched {Ak = 0) when the pump field is polarized along 
the slow axis and both signal and idler fields are polarized along the fast axis of the birefringent 
waveguide. 

The spectral entanglement generated by the Hamiltonian in Eq.lSlcan be viewed as a conse- 
quence of energy and momentum conservation. To quantify these correlations, we rewrite the 
Hamiltonian in terms of a minimal set of broadband modes via the Schmidt decomposition BOl 

N 

HsmM= Y, c,nAl{co,)Bl{cOi) + H.c. (4) 

where AJ,(®.v) = I dcOs^m{cOs)a\cOs) and B},^{cOi) = JdcOiy/i„{cOi)d''{cOi). The set of functions 
{£,,„}, {Wm} are called Schmidt modes, which define an orthonormal basis for the signal and 
idler Hilbert spaces respectively. This decomposition shows the evolution due to SFWM is 
equivalent to an ensemble of two-mode squeezing operators e'^sFWM=5^^^^ ^Sa2,B2'^--- where 
Sa„,,b„, is a two-mode squeezing operator on modes A„, and B^ ll2n . 

In general, the Schmidt decomposition in Eq.fflincludes significant contributions from mul- 
tiple modes so that A^ > 1. As previously discussed, a detector that cannot resolve these dif- 
ferent frequency modes leads to heralding of mixed quantum states. To restore high purity, 
one can employ spectral filters to remove higher modes but at the cost of reduced count rate. 
In contrast, the approach we adopt is to design the source to emit only in a single pair of 
Schmidt modes lfT6l [iTl IT2l 1411 . This allows heralding of high purity states without filtering. 
We use the singular value decomposition to numerically perform the Schmidt decomposition in 
Eq.H This allows one to predict the heralded photon purity given source parameters a{co) and 

To investigate the frequency-mode structure of our source, we measure the marginal photon 
number distribution. An ideal single-mode source would exhibit a thermal distribution, while 
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Fig. 3. Unheralded ^^ of the signal field shows control over the amount of spectral en- 
tanglement between the signal and idler fields as Kkp is adjusted (red). Theoretical curve 

(dashed) and data (points) are shown. Filtering out the peripheral lobes in the joint spec- 
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tral intensity with a 4.5 nm filter on the signal field is calculated to give gj, = 1.98 at 

AAp = 3 nm (blue). Insets: joint spectral intensity (JSI) measurements demonstrate spectral 



entanglement control. The FWHM of the pump envelope \a\ 
function |(^| (purple) accurately predict JSI orientation. 



(white) and phase-matching 



a highly multi-mode emitter gives Poissonian statistics ll42l . In our source, the transition from 
single-mode to increasingly multi-mode behavior can be readily adjusted via the pump band- 
width tsXp fTSl. To demonstrate this, we measure the autocorrelation ^.m- (0) as tskp is varied, 
as shown in Fig. BI For AAp = 3.1 nm we measure our optimal giv (0) — 1 .86 ± 0.02 which is 
close to the ideal thermal result e'^' (0) = 2. Only APD dark counts are subtracted from the data 
used to calculate gL (0) in Fig. pi One can relate the number of excited Schmidt modes to these 
statistics using ^.L (0) = 1 + l/L™ km|^ = 1 +^ where P is the heralded purity. Thus, we have 
demonstrated P = 0.86 without any filteiing of the signal and idler modes. To our knowledge, 
no previous on-chip source has simultaneously demonstrated purities and efficiencies as high 
as P = 0.86 and rjp = 80%. 

The near single-mode emission of our source is further supported by joint spectral intensity 
measurement. A spectrally uncorrelated state has a factorable joint spectral amplitude, and thus 
a factorable joint spectral intensity. In the middle inset of Fig. [3] we find that the major and 
minor axes of the central ellipse are parallel to the A,-,, axes, which shows a factorable intensity 
I{Xi,Xs) = Ii{Xi)Is{Xs). As AAp is adjusted away from this optimal value, the joint spectral 
intensities in the left and right insets of Fig.[3]indicate entanglement with an increasingly tilted 

(2) 

central lobe. These spectral correlations are in agreement with the measured decrease in gls {0). 
At the optimal AAp, the slight deviation from ideal thermal statistics, and corresponding 
reduction in heralded state purity, is principally due to peripheral lobes in the phasematching 
function. These arise from the hard-edge boundaries of the waveguide and are faintly observed 
in joint spectral intensity plots. A 4.5 nm filter on A, would suppress the small lobes, which 

(2) 

would yield gss {0) — 1.98, and corresponding P = 0.98, while still passing 90% of photons. 
Such a filter leaves rjp unchanged, as it would only be applied to the heralding signal photon. 

Our source offers large spectral tunability of the signal and idler fields. Fig. |4] shows that 
SFWM is phase-matched over a wide range of pump wavelengths. The inset figures illustrate 
that over this entire range AAp can be adjusted to achieve factorable output. For larger Ap the 
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Fig. 4. Theoretical phasematching curves sliow thie signal (blue) and idler (red) wavelengths 
that satisfy the condition M = for a range of pump wavelengths. Insets show predicted 
SFWM joint spectral intensities. A factorable output suitable for heralded production of 
pure photons can be produced over a wide spectral range, from visible to telecom wave- 
lengths, by adjusting Ay, and AAp (left to right: AAp = 3.1, 7.5, 20 nm). These calculations 
assume L = 4 cm and A« = 10^ . 



group velocities of the pump and idler photons become comparable, which results in a large 
spectral bandwidth for the idler. This relatively simple route to broadband emission could be 
useful, for example in quantum optical coherence tomography 1431 . In contrast, many applica- 
tions, such as quantum memories, instead require narrowband emission. A theory for cavity- 
enhanced SFWM has been developed that allows emission down to MHz bandwidths |44|. Cav- 
ities can be implemented using the refined Bragg grating technology available in silica 145,461 . 

5. Birefringence homogeneity 

Our analysis so far assumes the waveguide, and thus the wavevectors kp^i, are constant 
throughout the source. Imperfections in this regard can diminish the purity of heralded pho- 
tons. Due to our operation of a solid-clad guide far from its zero dispersion wavelength along 
with the excellent material uniformity of commercial silica I47,..15il . the dominant source of in- 
homogeneity is the birefringence. The spectral output is extremely sensitive to variation in the 
birefringence since the phase-matched wavelengths depend critically on this parameter lfT2]| . 

We consider two models of birefringence inhomogeneity and their effect on heralded state 
purity in Fig.|5] The effects of gradual changes during fabrication, for example variation in the 
writing laser power or local environment, are modeled as a birefringence that varies linearly 
along the length of the guide. Rapid fluctuations, on the other hand, are described as a random 
variation in the birefringence of a specified mean and standard deviation. For both models, the 
resulting phase-matching function is found by numerically integrating (Oy, O),) = /q e'"^*^'''c/z, 
where Afe(z) has an explicit spatial dependence. The corresponding joint spectral amplitude in 
Eq. [3] is then used to determine both the heralded purity via the Schmidt decomposition and 
the joint spectral intensities in the insets of Fig. 15] These simple models suggest our measured 
purity, to two standard deviations (gray band. Fig. l5]), coiTesponds to a birefringence inhomo- 
geneity of 5ma.r(An) < 3 • 10"^. 

Undesirable birefringence variations in sources can reduce the quality of quantum interfer- 
ence. Hong-Ou-Mandel interference of two single photons from identical sources produces a 
visibility equal to the photon purities. In the ideal case of 5 (An) = (left inset of Fig.B), we 
calculated a heralded purity of 0.98 when using a 4.5 nm filter on the heralding signal photon 
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Fig. 5. The heralded state purity resulting from two models of birefringence variation are 
shown. Gradual variation (solid line) is approximated by a linear dependence An{z) = 
Ano + S (An) x z/L. Rapid variation (dashed) is modelled by randomly fluctuating bire- 
fringence of mean (An) = Ang and standard deviation 5 (An). Random inhomogeneity in- 
creases the heralded purity for small 5 (An) due to suppression of the peripheral lobes in 
the joint spectral amplitude. All results assume L = 4 cm and Ano = 10^ . Insets show rep- 
resentative joint spectral intensities. The gray shaded region indicates the 95% confidence 
interval of our measured heralded purity. 



arm which still transmits 90% of photons. For an inhomogeneity 5(A«) = 3 • 10^^, the heralded 
purity, and thus interference visibility, remains high at 0.95 while the filter transmission is now 
84%. Accounting for source inhomogeneity, current fabrication methods thus appear sufficient 
to obtain high visibility interference from heralded sources. 



6. Conclusion 

We have demonstrated an on-chip source of heralded single photons that achieves extremely 
low loss (rjp = 80%) and high output state purity {P = 0.86) without any single-photon filter- 
ing. To our knowledge, no previous on-chip source has simultaneously demonstrated such low 
loss heralding of high purity states. Achieving this performance relied on spectrally factorable 
photon emission which was enabled by the mode control allowed by source integration. An 
estimate of birefringence inhomogeneity suggests current fabrication methods are sufficient for 
high visibility interference between multiple sources on the same chip. 

Our source meets several loss thresholds for quantum-enhanced applications. Interferometric 
phase estimation, using single-photon sources with rjp = 80%, can achieve a precision better 
than any classical probe field with the same number of photons fl). For linear optics quantum 

(2) 

computing, the high rjp and low gpi {0) demonstrated here enables entangling gates to violate 
Bell inequalities without postselection IH. Exploiting this performance in future work is facil- 
itated by the silica-chip architecture shared by our source and many recent integrated quantum 
optics experiments lH [H E] |32l [33 EH [35J . 

In the longer term, multiplexing of sources either spatially with a switching network B8ll49l 
or temporally with quantum memories ||50l may provide a route to construct many-photon 



quantum states |[8][l4l. Even with multiplexing, rjp and P still bound the composite source 
performance. Therefore, optimizing individual source performance remains critical. 

Our choice of SFWM in silica was motivated by the desire to minimize source loss. One 
direction for future work is to build similar SFWM sources at telecom wavelengths, where sil- 
ica loss is even lower. As we have shown, the spectral flexibility of SFWM allows factorable 
states to be generated at these wavelengths with proper adjustment of AXp . Quantum optics 
experiments in this spectral region, supplied by silica SFWM sources, could investigate a vari- 
ety of fundamental scientific questions that can feasibly be tested with larger numbers of single 
photons. 
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Note: After preparing our manuscript, we learned about related advances in heralded 
single photon sources by SPDC in a pp-KTP waveguide l[51]l . 



